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Abstract 

We calculate the Cosmic Microwave Background anisotropy bispectrum on large angular scales in the 
absence of primordial non-Gaussianities, assuming exact matter dominance and extending at second 
order the classic Sachs- Wolfe result 5T/T = <I>/3. The calculation is done in Poisson gauge. Besides 
intrinsic contributions calculated at last scattering, one must consider integrated effects. These 
are associated to lensing, and to the time dependence of the potentials (Rees-Sciama) and of the 
vector and tensor components of the metric generated at second order. The bispectrum is explicitly 
computed in the flat-sky approximation. It scales as Z~ 4 in the scale invariant limit and the shape 
dependence of its various contributions is represented in 3d plots. Although all the contributions to 
the bispectrum are parametrically of the same order, the full bispectrum is dominated by lensing. 
In the squeezed limit it corresponds to /^L al = —1/6 — cos(2#), where 6 is the angle between the 
short and the long modes; the angle dependent contribution comes from lensing. In the equilateral 
limit it corresponds to f^f — 3.13. 
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1 Introduction 



The linear approximation to cosmological perturbations has been so far sufficient and extremely 
fruitful, at least on large scales, before non-linearities induced by gravity become significant. How- 
ever, the accuracy of observations is now reaching a level such that all second-order effects, naively 
of magnitude ~ (10 -5 ) 2 , may become relevant. This is particularly important in the context of pri- 
mordial non-Gaussianities: second-order effects are in fact expected to give a signal of order /nl ~ 
few, which is not far from the present experimental limits [1, 2, 3]. A large amount of work has 
been done to study Cosmic Microwave Background (CMB) fluctuations beyond the linear approx- 
imation, in order to make predictions for the temperature bispectrum. As a complete calculation 
of the bispectrum is a daunting task, people concentrated on specific effects which are expected to 
dominate in particular limits. The bispectrum generated by the correlation between lensing and 
the Integrated Sachs- Wolfe (ISW) effect has been studied in [4, 5]. The one coming from lensing 
and the Sunyaev-Zel'dovich effect has been studied in [5]. In [6, 7, 8] the bispectrum generated by 
perturbations in the recombination history has been calculated. Refs [9, 10] (see also [II]) focused 
on very short angular scales where the signal is dominated by the non-linearity induced by dark 
matter clustering. A systematic control of all second-order effects in the Boltzmann equations is 
currently under study: see [12, 13] and references therein. 

In this paper we calculate the CMB bispectrum in the limit of large angles, i.e. on angular scales 
larger than the one subtended by the Hubble radius at recombination {0 > 1°); we do this assuming 
perfect matter dominance. Important, although, as we will see, partial results were obtained in this 
regime in [14, 15, 16]. 

Our calculation can be seen as the extension to second order of the classic Sachs- Wolfe formula 

[17] 

where <3? e is the Newtonian potential at recombination, which gives the large-angle prediction for 
the spectrum of the CMB fluctuations. As it is well known, this formula describes the angular 
variation of the temperature without considering the dynamics of the photon/baryon plasma, but 
only the gravitational redshift of photons from the last scattering to us. Therefore, it describes 
correctly the CMB anisotropics only in the limit where the scales considered are well out of the 
Hubble radius at recombination: the same restriction will apply to our calculation. The Sachs- Wolfe 
formula further assumes that decoupling took place when the universe was matter dominated - 
neglecting the transition between radiation and matter domination - and that the universe is still 
matter dominated nowadays, neglecting the present acceleration. At linear order this simplification 
is very convenient as the gravitational potential stays constant during matter dominance. At second- 
order the gravitational potential is no longer constant but the second-order metric during matter 
dominance is known [18] and can be written analytically as a function of the large-scale inflationary 
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perturbations [16, 19]. 

Clearly, these approximations do not hold in the real universe. However, our calculations give 
the exact bispectrum in the same limit in which the Sachs- Wolfe formula becomes exact: zero 
cosmological constant, recombination that happens much after equality and in the limit in which 
all scales are much larger than the horizon at recombination. This last limit can be imagined by 
thinking about an experimentalist making measurements in the far future, when the angle subtended 
by the Hubble radius at recombination is minuscule. The fact that our results become exact in a well 
defined physical limit is quite important, as on large angular scales the separation among different 
effects is in general gauge dependent. Therefore, one has to be careful in making approximations 
because neglecting some effects leads, in general, to a gauge dependent result. Besides its theoretical 
interest, we expect our result to represent a fair approximation to the real universe on large angular 
scales and it can be taken as a starting point for more elaborated calculations. 

Motivated by inflation, we assume that there are no vector or tensor perturbations in the initial 
conditions on super-Hubble scales. We perform the calculation of the CMB anisotropies by integrat- 
ing the photon geodesic equation during matter dominance using the so called "generalized Poisson 
gauge", which generalizes at second order the standard Newtonian gauge. Besides the Newtonian 
and curvature potential, at second order new terms are present in the metric, generated by the 
product of linear fluctuations: a vector mode in the dx l dt entry of the metric, and a tensor mode in 
the spatial part. 

All these terms contribute to the final CMB anisotropy. The time independent parts of the 
gravitational potentials give rise to second-order terms evaluated at last scattering, in analogy with 
eq. (1.1); their contribution was calculated in [14]. However, at second-order there are also terms 
integrated along the photon trajectory, similarly to what happens at first order when we depart 
from matter dominance with the ISW effect. The time-evolution at second order of the gravitational 
potential on sub-Hubble scales generates the well-known Rees-Sciama effect [20, 21]. But also the 
vector and tensor part of the metric contribute with two integrated terms. 1 All these terms contain 
a number of spatial gradients higher than the intrinsic terms, so that one may think that they 
can be neglected on large scales as suppressed by positive powers of k/(aH) at recombination. 
However, this conclusion is too hasty: these terms are integrated along the photon trajectory while 
modes progressively reenter the Hubble radius. Thus the ratio k/{aH) should not be evaluated at 
recombination but when the terms contribute to the time integral. We will see that all the integrated 
pieces give a contribution of the same order as the intrinsic terms in the equilateral limit. Actually 
the separation between intrinsic and integrated effects has no physical, gauge invariant meaning. 
For example, a part of the integrated vector contribution will turn out to be a boundary term. 

Another integrated contribution is gravitational lensing, due to the gravitational deflection of 

1 The integrated tensor contribution has been taken into account for the large scale anisotropies in [16]. 
The vector contribution has, to our knowledge, always been ignored. 
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the photon trajectory with respect to the line of observation. Although the effect of lensing on 
the bispectrum through its correlation with the ISW effect is well known [4, 5] (but absent in our 
calculation as we are assuming perfect matter dominance), we will see that lensing is important 
also when correlated with intrinsic contributions at last scattering. In particular, we will find that 
lensing gives a squeezed limit contribution of the same order as the one due to intrinsic effects, but 
which depends on the angle between the long and the short modes. The effect of lensing on the 
bispectrum was studied in [22] with the conclusion that its effect is suppressed in the squeezed limit 
by the tilt of the spectrum. We will see that this conclusion is not correct. 

In computing the CMB bispectrum we will employ the flat-sky approximation, which is valid 
for small angles of view. Given that at the same time we are interested in angles which are much 
larger than the Hubble radius at recombination, there is a quite narrow range of scales, 1 <C I <C 
Zistpeak) where our approximations hold. However, the flat-sky approximation greatly simplifies the 
algebra and makes the result much more transparent. The results will be given by 2-dimensional 
kernels B(l\,l2,h), which can be thought of as the 2d observable analogue of the kernels used (in 3 
dimensions) to describe the shape of the primordial non-Gaussianity [23]. 

The paper is organized as follows. In the next section we give the second-order metric in matter 
dominance in the generalized Poisson gauge as a function of the inflationary initial conditions and 
we calculate the temperature anisotropy at second order integrating the photon geodesic. In section 
3 we make a general discussion about the bispectrum of the temperature anisotropy in the flat-sky 
approximation and we calculate this quantity induced by a primordial non-Gaussianity of the local 
and equilateral kind. These are useful for comparison with our results. In section 4 we calculate the 
bispectrum using the results of section 2. The calculation is split (for convenience, not because the 
effects are physically distinguishable) in various pieces: intrinsic effects at last scattering, integrated 
vector contribution, integrated tensor contribution and lensing. The resulting total bispectrum is 
discussed in section 5 and conclusions are drawn in section 6. The flat-sky approximation is discussed 
in appendix A, while the details of the calculation of the Rees-Sciama effect are presented in appendix 
B. 

2 Second-order temperature anisotropics 

In this section we calculate the CMB temperature anisotropy at second order in perturbations, in 
the large angular scale limit and for matter dominance, as a function of the angle of observation. 
On large angular scales, the effect of second-order perturbations on the CMB fluctuations have been 
studied more generally in [24, 25]. Although we will later use the flat-sky approximation, the results 
of this section hold also in a full-sky treatment. 
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We are interested in the CMB temperature fluctuations, 

ST = T r M-% 

T T Q v ' 

where T {n) is the observed photon temperature in the angular direction n (n 2 = 1) and T Q is its 
average over the sky. For a black-body spectrum the observed temperature T G (n) is related to the one 
of emission T e (x e ) by Liouville's theorem: as phase space density is conserved in the propagation 
of photons (assuming there is no further scattering), the phase space density received in a given 
direction h is the same as the one at emission but with a temperature [26, 17] 

To (n) = ^T e (x e ) , (2.2) 

where uj e and lo q are the frequencies at emission and observation of a given photon. Notice that 
this statement is exact and therefore holds at any order in perturbation theory. In general, also 
the temperature at emission will not be isotropic, but will depend on the angle of emission. This 
dependence can be however neglected in our case, as we are interested in perturbations which are 
much longer than the horizon at recombination. 

We work in the so called generalized Poisson gauge and use conformal time r. In this gauge, the 
metric reads [18] 

ds 2 = a 2 (r) {-(1 + 2$)dr 2 + 2u i dx i dr + [(1 - 2¥)<S ij + 7 ij ]dxMx- ? '} , (2.3) 

where oj{ is transverse, uj^i = 0, and 7^ is transverse and traceless, 7^ = = 7^. In the mat- 
ter dominated era, assuming that the amount of primordial gravitational waves is negligible, the 
components of this metric are [18, 16, 19] 



$ =<j) + [(f) 2 + d- 2 {d j( t>) 2 - sartdidjidiWjt)] 
2 

21a 2 H 2 



+ d ~ 2 [md^f + 5(d 2 <P) 2 + 7d l( pd t d 2 4>] , (2.4) 



(,2 + V 2 (<^) 2 - 2d-%d j {d i ct>d j c 



* =(/)- 

+ l^TJp d ~ 2 [2( W) 2 + 5(d 2 <P) 2 + Id^d 2 ^] , (2.5) 

coi = - ^d- 2 [d 2 ^ - d^didjid 2 ^)] , (2.6) 

lij = - 20 Q - ^) d- 2 P^ kl (dM4>) . (2.7) 

The scalar quantities $ and ^ are the Newtonian and curvature potentials, respectively, while we 
will refer to Wj and 7^ as the vector and tensor components of the metric. The metric is expressed in 
terms of (j>, the time-independent quantity representing the initial curvature perturbation generated 
during inflation. Indeed, <p is simply proportional to the (non-linear) curvature perturbation on 
uniform density hypersurfaces (: on super-Hubble scales, where £ is constant, 

4> = -|c (* « o,H) . (2.8) 
5 



In the following we are going to assume that £ on large scales, and therefore (f>, is perfectly Gaussian, 
which is a very good approximation for example in minimal single field inflationary models [27, 28]. 
In the expression for tensor modes, the spherical Bessel function j\ (x) is given by j\ (x) = sin(x) /x 2 — 
cos(x)/x, while -Pj^ is a transverse traceless projector defined as 

PSl = \ (PikPjl + PjkPil - PijPkl) , (2.9) 
where Pjj is a symmetric transverse projector given by 

Pa - <% " 9 ~§- ■ (2-io) 

It can be expanded to give 

p?l {d k <f>di<f>) = -d- 2 [a 2 e <% + didjOo + 2(d 2 ( pd i d J 4> - d^d^)] , (2.11) 
0o = -^- 2 [(aW-(W) 2 ] . (2.12) 

In order to study the photon redshift we must solve the photon geodesic equation from last 
scattering to us, taking into account the perturbations of the metric above. The photon geodesic 
equation can be written as 



with 



dP M 



\d^P a P^ (2.13) 



dA 2 

where P^ = dx^ /d\ is the four-momentum of the photon, P^P^ = 0. The frequency of a photon with 
four-momentum P M as measured by an observer with four- velocity -u M , is given by u> = —P^u^. For 
simplicity, we choose the observer today to have zero spatial velocity, u % = 0. Indeed, any peculiar 
motion of the observer leads to a dipole anisotropy that can easily be subtracted. Furthermore, since 
we are interested in the large angular scales, we neglect also the Doppler effect due to the velocity 
of the photon/baryon fluid at recombination, which vanishes on super-Hubble scales at decoupling. 
Thus, we choose also the emitter to have zero spatial velocity, u\ = 0, so that we have u) = —Pqu° 
both for the observer and the emitter. Making use of the normalization condition of the four- velocity, 
u^u^ = — 1, one obtains uj = —Po/^—goo, and thus 

loo = Pq(t ) y-goole > 2 M , 

uj £ Po{r e ) ^-g o\o ' 

In order to compute Pq we need to solve the time component of eq. (2.13). Using that P° = dr/dA 
and plugging the metric (2.3) into eq. (2.13) yields 

P°^= ng a/3 P a pP - aV(P ) 2 + a^pVp* + a 2 ( - + \l[)P l P j , (2.15) 

where by a prime we denote the partial derivative with respect to conformal time, ' = d/dr, and 7i 
is the conformal Hubble rate, TC = a' /a. One can immediately notice that the first term on the right 
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hand side vanishes because of the massless condition P M P M = 0. Note also that, as we are studying 
perfect matter dominance, the two potentials <E> and ^ are constant at linear order, see eqs. (2.4) and 
(2.5). Thus, their time derivatives <&', ^>', together with uj\ and 7^, are all second-order quantities. 
One can therefore replace the zeroth-order expression P % = — P°n* into this equation; furthermore, 
using the background relation P° = —Po/a 2 , the geodesic equation can be finally rewritten as 





that upon integration yields 

Pq(To) 
Po(Te) ~ 



Plugging this expression into eq. (2.14), one obtains the photon redshift up to second-order as a 
function of the metric perturbations, 



lo a e / 1 + 2$ 



uj e a Q V 1 + 2$ 



l + ^ T °dr (V + VJr'+u;^-^ 



(2.18) 



Now we need to relate T e (x e ) on the right hand side of eq. (2.2) to the metric perturbations at 
decoupling. Since we concentrate on large angular scales, we only need the super-Hubble relation. 
We will use adiabatic initial conditions. In this case the dark matter energy density p rn simply scales 
as the third power of the temperature, 

Prn OC T e 3 . (2.19) 

In the matter dominated era, the energy density of dark matter is related to the metric perturbations 
through the Einstein equations, in particular through the energy constraint equation. On super- 
Hubble scales, i.e. neglecting spatial gradients, and using the fact that the potentials <!> and ^> at 
first order are time-independent in the matter dominated era this reads, up to second order, (see for 
instance eq. (196) of [29]) 

3H 2 (1 - 2$ e + 4$2) = 8vrGp m , (2.20) 

where H is the Hubble rate. Using the background Friedmann equation and eq. (2.19) above, this 
equation can be rewritten as 

T e = (l-2$ e + A<S> 2 e ) 1/3 f e , (2.21) 

where T e is the average temperature at emission, which simply scales as the inverse of the background 
scale factor. 

This equation can be derived in a simpler way [15] taking into account that, at recombination, all 
the modes that we are considering are much longer than the horizon and adiabatic. This means that 
each local observer will see a completely unperturbed history at any order in perturbations. Indeed, 
the vector and tensor components of the metric, eqs. (2.6) and (2.7), are suppressed by powers of 
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k/(aH) and can be neglected at recombination. 2 The same holds for the time dependent part of $ 
and i.e. the second lines of eq. (2.4) and (2.5). This means that the metric on large scales takes 
the form 

d s 2 = a 2 (r) {-(1 + 2$)dr 2 + (1 - 2* )S ij dx i dx j } , (k < aH) , (2.22) 

where <3? and ^ are now time independent and slowly varying in space. Locally, i.e. on scales of 
order of the horizon at recombination, this metric describes an unperturbed universe as the terms 
with $ and ^ can be taken to be constant in space and reabsorbed with a change of coordinates. 
In particular, the evolution is unperturbed in terms of a new conformal time f which satisfies (in 
matter dominance aar 2 ) 

r 4 (l + 2$)dr 2 = f 4 df 2 , (2.23) 

i.e. f = t(1 + 2<£) 1 / 6 . The temperature on a r = const, surface will be perturbed, because the same 
value of r corresponds to different moments along the unperturbed evolution, i.e. to different values 
of f. As T oc 1/f 2 we have 

T e = (1 + 2<S> e )~ 1 / 3 T e j ( 2 2 4) 

which coincides, at second order, with eq. (2.21). 

Now, let us plug both this equation and eq. (2.18) into eq. (2.2) and write the observed CMB 
temperature up to second-order as a function of the metric perturbations, 



T ' tF ' ) = t T '\lTTWj 1 + 2 ' l '' ) 



1 + rdr($' + *' + Jifj - 

J Te 



(2.25) 



Note that on the right hand side of this equation, the gravitational potential at the observer, ^*o? 
does not depend on the direction of observation. Thus, its dependence can be simply reabsorbed into 
the definition of T Q . Expanding this equation up to second order in the perturbation and plugging 
the right hand side in eq. (2.1) we finally obtain the CMB temperature anisotropics, 

ftT 1 5 l' r ° / 1 \ 

-^(n) = -d> e - -0> 2 + Jf dr (V + ¥' + - - V y nWj . (2.26) 

The first two terms on the right hand side of this equation have to be evaluated at the position 
of the emitted photon, x e . Since the second term is second-order, it can be simply evaluated at 
the background position nD e , with D e = t d — r e . Also the integral is second-order; thus it can be 
computed along the background photon trajectory, i.e. x(r) = nD(r), D{r) = t — t. However, the 
first term on the right hand side is a first-order quantity. Thus, at second order it must be evaluated 
at the perturbed position of the photon at emission. Expanding around the background position 
hD e we can write it as 

$(x e ) = <S>{hD e ) + 5x e ■ V<P{nD e ) , (2.27) 



2 Notice that in eq. (2.7) the prefactor in parentheses, 1/3 — ji(kr)/(kT), goes to zero for fcr — > 0, i.e. when 
the 7 mode is out of the horizon. 
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where 5x e = x e — nD e is the deviation from the background trajectory and we have used that $ = <fi 
at first order. 

In order to find 5x e we must solve the spatial component of the geodesic equation. Since V</> is 
already first-order we need to compute 5x e at first-order only. Thus, equation (2.13) gives 

= _2a 2 d^(P ) 2 , (2.28) 

dr 

where we have used that <I> + VP = 20 at first order. This equation can be integrated using the 
background relation P° oc 1/a 2 . The spatial gradient can be decomposed along and orthogonally to 
the background photon trajectory. Since <f) is time-independent, the component along the photon 
trajectory is a total derivative. Furthermore, raising the spatial index with the first order metric 
and then using P° oc (1 — 20) /a 2 one obtains 

— = — = -V(l + 2<P) + 2j dr'V\4>, (2.29) 

where we have defined Vj| = (J 1 - 7 — n l W)dj as the spatial gradient orthogonal to the line of sight 3 
and we have absorbed the dependence on <p a in the first-order definition of h, n % = — PI/ P® (1 + 20 o ). 
Integrating this equation and subtracting the background value hD e , after an integration by parts 
in the second integral one obtains the geodesic deviation 

?T 



5x, 



[To [To 

= 2n dr0-2 / 6Lt(t - T e )V \\<f> . (2.30) 



The first term on the right hand side, longitudinal to the line of sight, is the so-called Shapiro 
time-delay. This effect was discussed in [30] and we will discard it from the following discussion 
where we will concentrate on modes much shorter than the present Hubble radius, where the flat- 
sky approximation is valid. Indeed, since the integral of (f> tends to average to zero unless the mode 
wave- vector is orthogonal to the line of sight, it gives a negligible contribution to the CMB anisotropy 
for I 3> 1. The second term is the transverse deviation from the background trajectory, responsible 
for the lensing effect [31]. 

Including the lensing effect by re-expressing <3? e using eq. (2.27) and re-writing <E> in terms of 4> 
using the large-scale limit of eq. (2.4), i.e., $ = + 4> 2 + d~ 2 (dj<j)) 2 — 3d~ 4 didj(di<f)dj(/)), eq. (2.26) 
can be finally written as 



ST 



+ r dr (V + V + uft - ^tifj) + U ■ Vn4>e , (2.31) 



3 Notice that the direction perpendicular to the photon trajectory is parallel to the flat sky, so that, in our 
notation, the gradient is parallel to the sky. 
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where a is the deviation angle given by eq. (2.30) as 

a = -2 ! °dr r ~ Te V||0. (2.32) 

On the right hand side of eq. (2.31), the subscript "e" means at the background position of the 
emitted photon, hD e . The first line of eq. (2.31) was found in [14]. It represents an intrinsic effect 
due to the combination of the Doppler effect and the adiabatic temperature fluctuation of the plasma 
at recombination. The second line contains the Rees-Sciama effect, due to the second-order time 
evolution of the scalar potentials, and the effect of the time dependence of the vector and tensor 
components of the metric. Finally, the last term in the second line of eq. (2.31) represents the lensing 
effect. All these effects were discussed for a more general metric in [24, 25]. 

There is a nice way to check the factor (pg/18 in the expression (2.31) which, as we will see, is 
important for the squeezed limit of the bispectrum [14]. Let us take the limit in which one of the 
two Fourier modes of the initial conditions <fi e becomes infinitely long. This mode is still out of the 
horizon today and therefore cannot affect any physical observable. Let us check that this is indeed 
the case. When one of the wavevectors goes to zero, all the terms containing spatial derivatives in 
the expression above vanish, as it is clear from the explicit form of the metric eqs. (2.4)-(2.7). One is 
left only with the first two terms which, up to second order, it is useful to rewrite in an exponential 
form [15] as 

jct> r i 1 i 

~ e 0e/3 - 1 . (2.33) 



ST 



1 j. 1 
3 + l8 ( 



J e 



At first sight it looks as if the constant mode could affect observations through the second order 
term, which mixes a short mode with the constant one. This actually is not the case as the constant 
mode also affects the average measured temperature. Indeed the well defined measurable quantity 
is given by 

T (n) -% _ 

T " <e*-/3) • [2 ' 64) 
Now we see that indeed a constant contribution to <f> e cancels out: the quadratic term cancels with 
the redefinition of the average temperature. Notice that this is only possible because of the exact 
numerical coefficient 1/18 in front of the quadratic term. For the calculation of the bispectrum we 
are only interested in modes inside the Hubble radius at present time, thus it is not necessary to 
modify eq. (2.31) to take into account the correct average temperature as in eq. (2.34). 

In this way we also understand why the argument presented in [22] for the squeezed limit of the 
3-point function is not correct. In that reference it is argued that a term like (f>e/18, which induces 
a correlation between short and long modes, cannot exist, as it would imply - as in eq. (2.33) - that 
a mode which is still out of the horizon gives a measurable effect. What was neglected is that the 
same mode would change the average of the measured temperature. 
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3 The CMB bispectrum and its shape 

In this section we will discuss the CMB bispectrum and its shape dependence. We will use the 
flat-sky approximation. Even though this approximation is not very good for the lowest multipoles, 
the expressions that we will derive are much more transparent than using a full-sky treatment. 

In the flat-sky approximation (see appendix A) the Fourier transform in the sky of the temper- 
ature anisotropies is 

a r = J d 2 m^(n)e- r "\ (3.1) 

where we have decomposed h into a part orthogonal and parallel to the line of sight as n ~ (rh, 1) 
(see appendix A). The spectrum of the 2-point function is defined as 

{a t a F ) = (2ir) 2 6(T+?)C l . (3.2) 

We can rewrite the standard linear Sachs- Wolfe term in eq. (2.31) in Fourier space, 

£<*> = / ^5*^- (3 - 3) 

As explained more accurately in appendix A, it is convenient to separate k as the sum of a 2- 
dimensional vector parallel to the flat sky and a component orthogonal to it, 

k = (k\\k L ). (3.4) 

Using this decomposition and inserting eq. (3.3) in eq. (3.1) one obtains 

ar= J Wf \h^ De (2vr) 2 5(f- BD e ) . (3.5) 
From this expression the power spectrum defined in eq. (3.2) reads, 

a - 55* . < 3 - 6 > 

where for simplicity we have used a scale invariant power spectrum for the gravitational potential </>, 

<^^) = (27r) 3 <5(fc + fc')p. (3.7) 

We are interested in the ensemble average of the product of three ap Thus, we define the CMB 
bispectrum B(li,l2,fa) as 

{a f a m ) = (2w) 2 S{h + f 2 + l 3 )B{l 1 ,h, f 3 ) . (3.8) 

Translational and rotational invariance reduce the number of degrees of freedom of B to three 
independent variables only, for instance h,l2,h- This is completely general, but in the particular 
limit that we are studying (large scales and perfect matter dominance) we will also see that the 
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leading contributions to the bispectrum are scale invariant, i.e. the amount of non-Gaussianity is the 
same at long and short scales. Mathematically this implies that the function B is a homogeneous 
function of degree —4, 

B(Xl 1 ,Xl 2 , Xl 3 ) = \- 4 B(h,l 2 , f 3 ) , (3.9) 

which further reduces the number of degrees of freedom to two, for instance the ratios r 2 = h/h 
and r 3 = h/h- Without loss of generality we can assume < r 2 < r 3 < 1; the triangle inequality 
implies r 2 > 1 — r 3 . This is very similar to what happens when one studies the shape dependence of 
the primordial 3-point function of the curvature perturbation [23], with the difference that here we 
are in two and not three dimensions. 

We are interested in the dependence of B on the two ratios r 2 and r 3 , which describes how the 
bispectrum changes as we change the shape of the triangle in Fourier space. The possibility to mea- 
sure a bispectrum depends on its signal to noise ratio S/N, which is given in flat-sky approximation 
by [32] 

1 fd 2 l 2 d 2 l 3 B(h, l 2 ,l 3 ) 2 



The overall scaling in / is fixed by eq. (3.9) and (3.6): the integrand scales as l~ 2 . To study the 
shape dependence one can look at the quantity 

r 2 r 3 B(l,r 2 ,r 3 ) . (3.11) 

The square of this quantity is in fact proportional to the integrand in the expression above and thus 
quantifies the contribution to (S/N) 2 of triangles with a given shape. To be more precise one could 
rewrite the expression (3.10) for (S/N) 2 as an integral over the two ratios r 2 and r 3 

2 

(3.12) 



(S/N) 2 oc J dr 2 dr 3 



„3/2 3/2 

{ B(l,r 2 ,r 3 ) 



(2r 2 + 2r| + 2r 2 r 2 - 1 - rf - rf)V4 ' 



Therefore it would seem appropriate to consider the function in brackets as a measure of the S/N 
contribution; in this way in fact the integral of the square of the function over an r 2 , r 3 region would 
directly give the contribution of those shape configurations to (S/N) 2 . This would exactly parallel 
what is done in [23] to study the shape dependence of the primordial 3-point function. However in 
this way we would introduce a spurious divergence in the plots for flattened configurations when all 
the sides of the triangle are aligned: indeed, the denominator of the expression above blows up in 
this limit. This is just a consequence of describing the triangle shape in terms of r 2 and r 3 and it 
does not imply that flattened triangles are indeed more important. For this reason we prefer to plot 
f 2 r 3 B(l,r 2 ,r 3 ) in the following. 

For comparison with the results that we will derive later, it is interesting to study the function 
(3.11) when the CMB bispectrum is dominated by a primordial contribution. Two interesting cases 
are given by the so-called local and equilateral shapes [23] . 
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3.1 The local shape 

A popular shape, usually used in data analysis, is the one obtained when the potential <f> contains a 
non-linear correction in coordinate space, 

m = 4> 9 (S) - f^^x) - (4> 2 g )) . (3.13) 

(We are using the same sign convention for f^ 81 as Komatsu et al. [33].) In this case, the 3-point 
function of the gravitational potential <f> is 

<^AAs> = (2-) 3 ^(fci + k + k 3 )(-2f^A 2 ) (-I3 + -I3 + -I3) . (3.14) 

If the non-linear correction (3.13) dominates over those computed in the previous section, then 
can be simply computed using eq. (3.5). By taking the ensemble average of the product of three aj- 
and using eq. (3.14), the CMB bispectrum induced by local non-linear corrections reads 

? local _ ^1,1,1 



B ° c& - -^-\M + m + w' (3,15) 

Note that by rescaling I2 and I3 we can pull out an overall factor /j~ 4 and rewrite this bispectrum in 
terms of the two independent variables r 2 and r^, 

R local _ 2 /NL a1 ^ 2 ( 1 , 1 1 \_\ filfCt 

~ 27--'/ \~ ' ~ ' TTTTI ) • (■> A( >) 



n r 2 r 3 



In the following we will always use this trick and plot the bispectrum as a function of r 2 and r% 
setting h = 1 and A = 1. The shape corresponding to eq. (3.16) is plotted in figure 1. 



3.2 The equilateral shape 

Another theoretically motivated shape for the primordial 3-point function is the so-called equilateral 
shape, that can be described by [23] 

%A-Aa> = (^te + & + ^)(-6/X^ 2 ) (-^13 " T^T2 + T^U + 5 P™) . 

(3-17) 

Notice that the divergence in the squeezed limit is in this case milder than for the local shape, due 
to a cancellation among the various terms. We can compute the CMB bispectrum similarly to what 
is done in the local case. It is convenient to define 

Vl = ki{D e /h) , y 2 = ki{D e /h) . (3.18) 
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Figure 1: The CMB bispectrum on large angular scales induced by primordial non-Gaussianities 
of the local form for /jJl 3,1 = 1- According to its definition, the bispectrum is negative for positive 
/nl^j thus, we have plotted it with an overall minus sign. 

With such a definition, using eq. (3.5) for the and eq. (3.17) for the expectation value of three 
gravitational potentials, one finally obtains 



ofcquil ,2 r +ca 

* equil = ^TW" / d Vl dy 2 
9(2vr)^ y^oo 



1 1 

+ 



2(yl + rl)3/2( y 2 + r 2 ) 3/2 3 (y 2 +r 2 )(y 2 +r 2 )((yi+2/2) 2 + 

1 

+ 5 perms. 



(yj + r\yi\yl + r 2 2 )(( yi + y 2 f + r 2 ) 3/2 

(3.19) 

Here and in the following we sum over all permutations of (r\, r 2 , r$) and we subsequently set r\ = 1. 
The integrals cannot be done analytically but the result is plotted in figure 2. 

From figures 1 and 2 we see that the CMB bispectra preserve in 2d the qualitative features of 
the primordial 3-point functions: the signal is peaked on squeezed and equilateral configurations 
respectively. 



4 Computing the CMB bispectrum 

In this section we compute the CMB bispectra due to the different second-order contributions in 
eq. (2.31). For comparison, we will use the two typical primordial shapes, local and equilateral, 
discussed above. We are only interested in computing the CMB non-Gaussianities generated in 
the Sachs- Wolfe limit; thus, as already mentioned, we will assume that there is no primordial non- 
Gaussianity, i.e. that the curvature perturbation on uniform-density hyper surf aces, is Gaussian 
on super-Hubble scales. Consequently, from eq. (2.8) it follows that 4> is Gaussian. 
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T'l 

Figure 2: The CMB bispectrum on large angular scales induced by primordial non-Gaussianities 
of the equilateral form for f^f^ = 1. According to its definition, the bispectrum is negative for 
positive /nl" 11 ) thus, we have plotted it with an overall minus sign. 

4.1 Intrinsic contributions at last scattering 

Let us start by computing the CMB non-Gaussianity due to the second-order effects in the first line 
of eq. (2.31), i.e., 



4^ + V 2 (O%0) 2 " Sd-Zdid^drf)) 

io o 



(4.1) 



This contribution has been first derived in [14] and its bispectrum and detectability have been 
studied in [34]. Note that, although we have dubbed it "intrinsic", this contribution is not physically 
separable from the other second-order contributions integrated along the photon path that we will 
study below. 

The momentum-independent quadratic term, gives a contribution to the bispectrum 

exactly of the local shape, equivalent to /nl 3 " 1 = —1/6 [14], in eq. (3.16). Its contribution does not 
vanish in the equilateral limit. We can compare it to an equilateral contribution by evaluating its 
bispectrum in the equilateral configuration. We find 

B-V 6 (1A,1) 



-0.24 , (4.2) 



where we have evaluated i? eqml (l, 1,1) for /nl" 1 = 1- We conclude that this contribution is equivalent 
to /nl"' — —0.24 in the equilateral limit. 

In order to compute the contribution from the momentum-dependent term, we rewrite it as 

1 ^)- 2 <F) ; rhA 2 - !) '!) ;) i!) r, ;) , , I = I d Pl d P2 f intr C^ rh- rhr, J@l+P2>* D e 



-d-\d i( j> e y - d-^djid^ed^e) = J ( "; ( "3 / lntr (pi,P2) fafaJ&W" , (4.3) 
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where / intr (pi,P2) is a kernel defined as 



1 Pi • P2 _ Pi Pi + (pj + P 2 )(Pl -P2 ) + (Pi ■ P2) 

3(p 1 +p 2 ) 2 {P1+P2Y 



Note that this kernel vanishes in the limit of either p\ or p 2 going to zero. Thus, we expect this 
contribution to be suppressed with respect to the local shape in the squeezed limit. 
The Fourier transform in the sky of this contribution is 

<*,-= / ^^f^^A)^<^^^ )D 'W 2 S{t- (4+4)D e ) ■ (4.5) 

To compute the bispectrum we can contract this contribution, which is quadratic in (f), with the 
product of two linear Sachs- Wolfe effects, whose are given by eq. (3.5). By doing so, evaluating 
the 4-point function of <\> using Wick's theorem and the definition of the power spectrum, eq. (3.7), 
summing over all permutations, and using the definition of the bispectrum, eq. (3.8), one obtains 



942 /-+00 



9(2vr) 2 Z 



( 23/12/2 + r 



(y? + rlf'^vl + r 2 2 )3/2 ^ 6 ( (yi + y2 )2 + r: 



3) 



4(y? + r 2 )(y| + r 2 ) + 2{y\ +r\ + yl + r%)(2 yi y 2 + r 2 - r\ - r\) + (2 yi y 2 + r| - r 2 - r. 



2\2 



4((yi + y 2 ) 2 + r 2 ) 2 



+ 2 cyclic 



(4.6) 



The integrals in the expression above can be integrated numerically. The final result for the bis- 
pectrum coming from this contribution is plotted in figure 3. Its contribution is equivalent to 

/NT 1 ^ 1-21. 

Notice that this bispectrum is suppressed in the squeezed limit with respect to the local case in 
figure 1. This, as discussed, is a consequence of the derivatives in eq. (4.1). Notice also that the 
suppression, in the limit r 2 — ► 0, is linear in r 2 as there is one derivative acting on each <j> in eq. (4.1). 
Thus, in the plots (which include a measure r 2 r^) the function goes to a constant. This constant 
depends on the orientation between the long wavelength mode and the short ones as it is clear from 
eq. (4.4): indeed, in the figure we see that the limit r 2 — > depends on the direction from which the 
limit is approached. Notice that this behaviour is different from the case of primordial equilateral 
non-Gaussianity where there is a suppression going like r 2 in the squeezed limit - B eqml diverges 
logarithmically for r 2 — > 0, see eq. (3.19) - so that the plot in figure 2 goes to zero. Indeed, in this 
limit the 3d kernel (3.17) is suppressed by two powers of £3 with respect to the local shape and this 
behavior is typical of all equilateral models [35, 23]. 



4.2 Contribution from the Rees-Sciama effect 

At second-order in the perturbations, the Newtonian and curvature potentials $ and ^ have a 
constant and a time-dependent part. While the constant part given in the first line of eqs. (2.4) and 
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Figure 3: The CMB bispectrum induced by the momentum dependent intrinsic contribution in 
eq. (4.1). 



(2.5) dominates on large scales, on sub-Hubble scales one recovers the standard Newtonian limit 
[36], i.e. the two potentials become equal, $ = ^, and grow as the scale factor, $ oc (aH)~ 2 oc a, 
where we have used a oc r 2 . Thus, we expect the photon frequency to be affected by an integrated 
effect. This is the so-called Rees-Sciama effect [37], given by 



ST, 



n 



(4.7) 



Its contribution to the CMB bispectrum has already been considered in [20, 21] although these 
analysis were restricted only to the diagonal terms of the bispectrum. More generally, the bispectrum 
from the Rees-Sciama effect has been studied in [38]. 4 

Symmetrizing over the momenta, we can rewrite the integrand in eq. (4.7) as 



+ = J_ / d3 Pi d ^ 2 
D e J (27r)3(2vr)3 



/ RS (P1,P2)' 



'Pi Yp 2 1 



J(pi+P2)-nD(r) 



(4.8) 



where / RS is an explicitly time-dependent kernel derived from eqs. (2.4) and (2.5) defined as 

4(pi • V2? + Wpjpl + 7(pj +pl)(pi ■ P2) 



f RS (Pl,P2) 



-tD p 



(4.9) 



Note that we have multiplied it by D e to make it dimensionless. The Fourier transform in the sky 
of this contribution is given by 



dr C d?pi d 3 p 2 
WJ (2vr)3 (2vr)3 



/ RS (A,P 2 ) ^ 1 ^- 2 e i( ^ + ^ )D(T) (2vr) 2 ( 5(r- (p{ +p\)D{r)) . (4.10) 



'As there is an error in the derivation of eq. (23) of [38], our results cannot be compared with that reference. 
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As done for the intrinsic contribution, in order to compute the bispectrum we need to contract 
dj- in the above equation with the product of two linear Sachs- Wolfe effects, whose are given 
by eq. (3.5). Note, however, that the Rees-Sciama kernel / RS in eq. (4.9) is higher order in the 
spatial gradients with respect to the intrinsic kernel / mtr of eq. (4.4), so that one may think that its 
contribution to the bispectrum will be relevant only on short scales. Indeed, since we are correlating 
the Rees-Sciama effect with the linear Sachs- Wolfe effect, which takes place at the last scattering 
surface, one may naively conclude that its contribution to the bispectrum is suppressed in the limit 
of large angles, i.e. in the limit where gradients are much smaller than the Hubble rate at decoupling. 
However, this is not the case. Indeed, the correlation with what happens at the last scattering surface 
does not vanish immediately for r > r e , but for a given mode /, it remains large for r < = D e /l 
and after that decays exponentially. In other words the correlation decays when the distance from 
the last scattering surface is of the order of the typical wavelength. In appendix A we explain better 
this point with a simple example. Now, since the Rees-Sciama effect grows with r, the contribution 
to the bispectrum will be maximal for r « t*. Using that k ~ l/D e one has that the maximal 
contribution comes for /ct* ~ 1 so that the gradients are not suppressed at r* and one expects the 
Rees-Sciama contribution to the bispectrum to be of the same order as one of the intrinsic terms. 
Notice also that, as for the intrinsic kernel, also the kernel (4.9) vanishes in the limit of either pi or 
P2 going to zero; we thus expect the Rees-Sciama bispectrum to be suppressed in the squeezed limit 
with respect to the local shape. 

Let us move to the explicit calculation. It is convenient to define 

X= (T-T e )(h/D e ) . (4.11) 

By contracting given by eq. (4.10) with the product of two linear contributions given by eq. (3.5), 
using Wick's theorem and the definition of the power spectrum, eq. (3.7), to rewrite the 4-point 
function of (f>, and summing over all permutations one obtains, by using the variables y\ and y2, 
2A 2 1 f' 1 r+oo r i 

/3 o 3 9 9 5. 9 9 . 5 (y? - yl + r? - rk) 2 \ 

V 2 2 2 2 (yi+y2) 2 + ri ) 

(4.12) 

Actually the result of the calculation is not proportional to (2ir) 2 5(li + h + h) as in the definition of 
eq. (3.8), but to (2tt) 2 5 + h) ^-r + ^) ano - permutations, as a consequence of the fact that we 
are correlating effects at different conformal times r. This is a bit surprising as the delta function 
8{h + h + h) is just a consequence of translational invariance. However, the discussion above implies 
that the bispectrum is exponentially suppressed when the triangle in Fourier space does not close, 
i.e. when h ~ 1. This can be checked explicitly in the expression (4.12). In appendix A we 
discuss a simple example in which this issue is made more transparent. 
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The above integrals are particularly challenging even numerically. However, some simplifications 
can be made. Since the integrand is exponentially suppressed for i > 1 by the rapid oscillations of 
e H2/i +U2)x ( one can pugh the upper limit of the integral in x to oo. Another simplification consists in 
neglecting T e {l\/D e ) in the first integral of eq. (4.12), which is justified by the fact that we consider 
only modes well outside the Hubble radius at recombination and thus T e {l\/D e ) ~ r e k <C 1. With 
these approximations eq. (4.12) can be rewritten as 



B 



RS 



2A 2 



1 



189(2tt) 2 If 



dx x 
5 



dyidy 2 e 



i(yi +V2)x 



1 



-t{ + -r l 2 + r£ + 2y m + -(y( + y z 2 



5(y? 



vl + rl 



(y2 +r 2)3/2( y 2 +r 2 )3 /2 

4? 



2 (yi + y 2 ) 2 + r\ 



2 cyclic 



(4.13) 



2 x 2 ^ ° 2 

We see that the bispectrum induced by the Rees-Sciama effect goes as I and it is parametrically 
similar to the intrinsic contribution discussed in the previous section. The analytical and numerical 
study of this expression is postponed to appendix B. The final result for the bispectrum is given in 
figure 4. 

As for the intrinsic contribution (4.6), in the squeezed limit r 2 — > the Rees-Sciama bispectrum 
is suppressed when compared with the local shape by r 2 , with a coefficient which depends on the 
angle. We show this analytically in appendix B. By comparing the Rees-Sciama bispectrum to the 
equilateral contribution, as we did for the intrinsic one, we find that the Rees-Sciama contribution 
is equivalent to f^f 1 — 0.74. 




Figure 4: The CMB bispectrum induced by the Rees-Sciama effect, eq. (4.7). 



4.3 Integrated vector contribution 

At second order, the non-diagonal part of the metric goi = a 2 Ui becomes non-vanishing and time 
dependent on sub-Hubble scales. Similarly to the time-dependent part of the gravitational potentials, 
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it induces an integrated effect on the photon redshift, given in eq. (2.31) by 

ST, 



-(h) D 



dr torn 1 . 



(4.14) 



As LOi is transverse we refer to this effect as the integrated vector contribution. As for the Rees- 
Sciama, to compute the bispectrum we need to correlate this integrated effect with the intrinsic 
temperature fluctuation at last scattering. Even though this effect is suppressed at last scattering, 
when modes are still out of the Hubble radius, it will give us a contribution to /nl of order unity, 
similarly to what happens for the Rees-Sciama effect. 
From eq. (2.6) we can rewrite the integrand as 

d 3 pi d 3 p 2 , V /- - 



u'h 1 = 



1 



D e J (2vr)3(27r) 



'Pl Yp 2 c 



i{pi+P2)-nD{r) 



(4.15) 



where / v is a kernel defined as 
2iD, 



f Y {Pl,P2) 



p\{n-p 2 ) + p|(n-pi) . _ ,2p\p\ + (pl +pl)(pi -P2) 
n ■ (pi + p2) — i 



(4.16) 



(pi + P2) 2 " ' ' (pi +P2) 4 

Note that the second term in the kernel (4.16) is proportional to h ■ (pi +f>2)- Thus, it is a time 
total derivative which can be trivially integrated in r in eq. (4.14). Therefore we have another term 
evaluated at last scattering, analogous to the intrinsic contributions studied in section 4.1, of the 
form 

XT 1 A 



^(n) D I [d-%{d* 



T 

This shows clearly that there is nothing really intrinsic about the contributions discussed in sec- 
tion 4.1: the splitting among the various effects is gauge dependent and only the total sum has a 
well defined gauge invariant meaning. 

One can then split the rest of the kernel orthogonally to and along the line of sight. Indeed, 
decomposing h into the parts orthogonal and parallel to the line of sight as h = (m, 1), the first 
term in eq. (4.16) can be rewritten as 



2iD P 



rh • {pjpl+pjpl) 
(pi + P2) 2 



+ 



PlP2" + P2P1" 

(pi +P2) 2 



(4.18) 



The first term of this expression is proportional to fa. Thus, it is higher order in l/l with respect to 
the second term and therefore negligible in the flat-sky approximation. Thus, the Fourier transform 
on the sky of the contribution (4.18) can be approximated with 

jj. 1 2 



2iD e 



ar = —- 



T ° <fr_ 
. ~D~ P 



I dp^dp2_ Plpt+Plpt , , ^pt + pi)D(r) (27r) 2 d( f_ ,J J )D(t)) 
J (27T) 3 (27T)3 {fr+fr? ^ ' { Wl+P^W)- 



(4.19) 
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Proceeding as in the case of the intrinsic and Rees-Sciama contributions, the total contribution 
from vectors can be written, using the variables y\ and y 2 , as 



B 



v 



4A 2 



(y 2 + r 2 )(2y 2 ( yi +y 2 



(yf+rf)3/2(y| + r |)3/2 



4 - 4)) + ivl + rl){2yi{ yi + y 2 ) - (r\ 



{)) 



/o 



2((yi + y 2 ) 2 + r|) 

)x (vi + r i)i/2 + (yj + 

2 
3 



(2/1 + y 2 ) 2 + 



+ 2 cyclic 



(4.20) 



The first piece, which is not integrated in x, comes from eq. (4.17), while the other term describes 
the contribution integrated along the line of sight. The integral over time can be dealt with as in 
the Rees-Sciama case: see appendix B. The final result for this bispectrum is given in figure 5. 
Again, the result is suppressed with respect to the local shape in the squeezed limit because the 
kernel (4.16) vanishes when either p\ or p 2 go to zero. The behaviour in this limit is qualitatively 
the same as in the Rees-Sciama case. This vector contribution is equivalent to /j^l"' — — 0.84 in the 
equilateral configuration. 



r 2 r 3 B^ 




Figure 5: The CMB bispectrum induced by the vector contribution, eq. (4.14). 



4.4 Integrated tensor contribution 

At second order, even in the absence of primordial gravitational waves, the part of the spatial metric 
not proportional to the identity matrix, a 2 jij, is non- vanishing and time dependent on sub-Hubble 
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scales. Thus, it induces an integrated effect given by 



— (n)D-y dr-^-nV, 



(4.21) 



which we expect to contribute to the bispectrum similarly to what happens for the vectors. As 7^ is 
transverse and traceless, we refer to this effect as the tensor contribution. From eq. (2.7) and using 
eqs. (2.11) and (2.12) to rewrite the transverse traceless projector, the integrand is 

_1 , = J_ f d 3 pi d 3 p 2 fTf-t -. 
2 % 3 nn D j (27r)3(2vr)3 

where the kernel / T is defined as 

5D, 



f T (Pl,P2)< 



PI Vf>2 1 



i(pi+P2)-h D(t) 



(4.22) 



/ (pi, £2) = -h{\Pi+P2\r)- 



(Pl -P2) 2 ~p\p\ 
(P1+P2) 4 



1 + 



(n • (pi +P2)f 
(pi +P2) 2 



2p{(n ■ p 2 ) 2 + 2p 2 2 (h ■ pxY - 4(pi • p 2 )(n ■ pi)(n ■ P2) 



(4.23) 



(P1+P2) 4 

and j'2 is a spherical Bessel function that appears from taking the time derivative of 7^, 

' h(kr)\' _ J2(fer) 
kr 



(4.24) 



As we did for the vector kernel / v , / T can be decomposed into a part parallel and orthogonal 
to the sky. The parallel part is higher order in 1/7 and thus negligible in the flat-sky approximation. 
Thus, the kernel can be approximated as 

{P1+P2) 2 * 



f (pi,P2) ^ — J2(|pi +P2\r, 



5D P 



(pi ■ P2) 2 - p\p\ 
(pi+P 2 ) 4 



1 + 



(pi +P2Y 
2p?(p 2 L ) 2 + 1p\(Pi) 2 ~ 4(pi • P2)pf L P 2 L 



(4.25) 



(P1+P2) 4 

The Fourier transform on the sky of this contribution is given by 

a ?= £ W e I ^^^^ hAp^ + ^ mT \^?5{l- (Pi +A)D(r)) . (4.26) 

With this simplification the time integral can be analytically computed and yields, expressing it in 
terms of the variables x, y\ and ?/2, 



° dx h(V(yi + y2) 2 + rix) ci(yi+y2)x = (2rf - (m + y 2 ) 2 ) 
x 6((yi +y2) 2 + rl) 



(yi + 2/2)73 coth ( yi+y2 — - J 



2((yi+y 2 ) 2 + r 3 2 )3/ 2 



(4.27) 
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plus an imaginary term odd under (2/1,2/2) — * {— Vi, —2/2) which does not contribute to the integral. 



This gives for the bispectrum 
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6((2/i + 2/2) 2 + 4) 



x 



X 



(2/f + rl)^(yl + r 2 )3/2(( yi + y 2 )2 + r |)2 

- 2 {ylivl + r 2 ) + y\{y\ + r 2 ) - 2/12/2(22/12/2 + rf - r 2 - r 2 )) 



(4.28) 



+ 2 cyclic 



and the final result is plotted in figure 6. Again, given that the kernel (4.25) goes to zero when 
either p\ or P2 go to zero, this shape is suppressed with respect to the local one in the squeezed 
limit. From figure 6 we see that the integrated tensor contribution is qualitatively similar to the 
intrinsic kernel, Rees-Sciama and vector contributions discussed previously. This contribution is 
equivalent to /nl"' — —0.61 for an equilateral configuration. 
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Figure 6: The CMB bispectrum induced by the tensor contribution, eq. (4.21). 



4.5 Lensing 

The deflection angle of a light ray as it propagates from the last scattering surface to us is given 
by eq. (2.32) (for a review of lensing effects on the CMB see [39]). For convenience we reproduce it 
here, 

a = -2 dr T ~ Te VU. (4.29) 
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The geometrical weight (r — r e )/(r — r e ) tells us that the effect is suppressed close to the last 
scattering surface. For this reason, usually the main contribution to the 3-point function due to 
lensing comes from the correlation of the photon deflection with the ISW [4, 5]. This effect is absent 
in our case as we are studying a universe with only matter. However, there is still the correlation of 
the intrinsic temperature fluctuation at last scattering with the lensing contribution given by 

ST 1 

— (n)D-a-V^e- (4.30) 

Similarly to the other integrated effects also this will give an effective /nl ~ 1- 

Let us compute the contribution to the bispectrum. Inserting the deviation angle (4.29) into 
eq. (4.30) and using = D e V||, the lensing contribution to the temperature fluctuation can be 
written as 

f ^ ^ i r dr (r - re) / • A)h^- hD{T) ^ D « ■ (4.3i) 

Taking the Fourier transform on the sky yields 

a r = 2 -£°dr(r-r e ) J ^3 •p1)0p 1 0p- 2 e^ D ^e^^(2^) 2 ,5(Z < - P 1^(r) -^D e ) . (4.32) 

As usual, we can compute the bispectrum by correlating this effect with the intrinsic temperature 
at last scattering. By doing so, we obtain 

f-oo roo ^,2 „2 

27(27 

Note that here one must sum over all permutations of (?"i, r 2 , r 3 ), including the anticyclic ones. The 
integrals above can be computed analytically, yielding 



Al roo roo '2 _ 'A _ 'l 

# lens = - — . 2/4 / dxx d yi dy 2 3 3 ^ 2 + 5 perms. . (4.33) 



8,4 2 



„2 Jl Jl 



This result is plotted in figure 7. Alternatively, this equation can be written as 

B '""=^ii +bpams - (4 - 35) 

Another method to derive the lensing CMB bispectrum is through the lensing potential ip defined 
as (see for example [39]) 

^(n) = -2 Pdr _ T ~J £ _ , <l>(h(T -T),T) . (4.36) 

Jre \ T o T e ){T T) 

The deflection angle (4.29) is obtained by taking the flat-sky gradient of this expression a = V^V- 
The correlation between the temperature at the last scattering surface and the lensing potential is 
given by 5 

( ^ } = ~h^r 2 1 dT Kl{h{T " Te)/De) 6{h + h{jo " T),De) ■ (437) 



5 As explained in [39], the divergence of the lensing potential at r D affects only the monopole, which can 
always be subtracted. 
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r 2 

Figure 7: The CMB bispectrum induced by the lensing contribution, eq. (4.30). 

The temperature fluctuation is localized at r e while the lensing becomes more and more important 
at later times. It is easy to see that the correlation is maximal at t* ~ t q /1, similarly to what was 
discussed for all integrated effects in section 4.2 (see also appendix A). 

As the integral is dominated by r <C to, we can approximate the 5 function with 5(l\ + 1%) to get 

(rP ti a r2 ) = (2v)H(h + l 2 )Cj^ , Cf^ = - i^. (4.38) 

The bispectrum can be written as [40, 32] 

£lens = _f x . f 2 ^ C T^ + ^ C T,^ + 2 permg > (4 g g) 

which coincides with eq. (4.35). 

Note that the expression in eq. (4.35) diverges in the squeezed limit. However, the form of the 
divergence depends on the direction one approaches the limit. One can compare the expression 
resulting from taking ri — > with the local form. This gives a contribution equivalent to /nl^ = 
— cos(2 6) where is the angle between I2 and l\ when one takes the limit. In the equilateral 
configuration, the lensing gives a sizable contribution, equivalent to — 2.87. 

It is possible to recover the lensing 3-point function in the squeezed limit in another way, which 
is physically more transparent and can be easily generalized to the case when the short wavelength 
modes are inside the horizon at recombination. We are going to calculate the 3-point function by 
first taking the long wavelength mode fixed and then studying its lensing effect on the short scale 
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2-point function. 6 At the end we average over the long wavelength mode. 7 Consider the 2-point 
correlation function of the temperature fluctuations in two different directions h\ and n 2 . In the 
presence of a long wavelength mode the real space 2-point function is lensed 

ST ST \ /5TST\^ _ , A _ A , . 

— (ni) — (n2)y = — — y [ni + a(ni)-n 2 -a(n 2 )J , (4.40) 

where we used the fact that the unlensed 2-point function just depends on the distance between the 
points. Obviously there is no effect if the two lensing angles are the same: the 2-point function is 
just translated. Expanding at first order and defining by fh\ and m 2 the components of hi and w 2 
parallel to the (flat) sky we have 



1ST ^ JT^ \ /STST\.^ _ . „ / ST ST \ 
^ — (ni) — (n 2 ) J = ( yf~f / i m i- m 2j + V^ — — )VjQ!i 



mi + m 2 



(mi -m 2 )j . (4.41) 



By assumption the lensing wave is of long wavelength so that we can evaluate the gradient of the 
lensing angle at the midpoint (mi + m 2 )/2. If we call fa = rh\ — m 2 , we have 



ST . , ST . \ /STST\ r . d /STST\ r . m 7 - m, 

— (ni) — (n 2 ) > = ( — — > m + — ( — — > m —±—± VjCXi 

T T '/tens \ T T / 1 dlogm \ T T / mm 3 



mi + m 2 



(4.42) 



We can now Fourier transform to Zi and / 2 . The result can be expressed in terms of Is = (h — Z 2 )/2 
and II = h + h-, where l and s stand for long and short wavelength, 

K^ie- = C ls + UlMIl) [ d 2 m-^— ( 5 -7^r%-\ H — — . (4.43) 

41 ' 2 s 7 dlogm \ T T / mm 

The 3-point function is obtained multiplying the above expression by ST jT of the long wavelength 
mode and averaging, 

(a ti a r2 a t3 } = (2,fS{h + T 2 + T 3 ) ■il Lj ( S ^a i y (l L ) [ d*m-±- [m] . (4.44) 

11 12 13 1 J dlogm \ 1 1 I mm 

The prime in the correlation between lensing and the temperature means that we have to remove 
the momentum conservation factor (2it) 2 S, which has been factored out. 

Let us evaluate the integral over rh, which describes the effect of lensing on the 2-point function. 
One may naively think that for a scale invariant 2-point function, which is the case that we are 
studying in this paper, the effect of lensing vanishes. Indeed, the calculations above are very similar 
to the ones leading to the consistency relation for the squeezed limit of the primordial 3-point 



6 It is easy to argue that the leading contribution in the squeezed limit is obtained when the lensing mode is 
of long wavelength. Indeed, lensing is effective far from the last scattering surface, but as we get far from it the 
correlation with the temperature fluctuation rapidly decreases. The loss of correlation happens at <~ t q /1, 

i.e. it is faster at high I, that is why the squeezed limit is dominated by a long lensing wave. 

7 This discussion is inspired by the derivation of the consistency relation for the squeezed limit of the 
primordial 3-point function in single field inflation [27, 41, 42]. In particular we will parallel the explicit 
derivation done in sec. 2 of [42]. 
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function [27, 41, 42]. In that case, however, the integral over fn does not contain the angular weight 
m?m % /m? . Without this terms the integral vanishes for a scale invariant spectrum: indeed the 2- 
point function in real space is a logarithm of the distance, so that its log-derivative is a constant. 
The Fourier transform of a constant is 5 {Is) which vanishes for any non-zero Is- 

The situation is different in the presence of the angular weight m?m % /m 2 . To be more explicit, 
let us introduce a scale dependence in the 2-point function and evaluate the integral in eq. (4.44) 
for a power spectrum of the form C\ = C • ;- 2 +("s- 1 ) j which corresponds to a 2-point function going 
as m~( ris ~ l \ to see that the result does not vanish for n s — > 1. The integral can be written as 



/ 



d 2 m-—- / — — \[ m ] — ^ e -** = -(„, - l)%gi / d 2 m e"^ 
dlogm \ T T / L J m m V V 2 J \T T 1 



= -(n s -l)4#CW~ 2+(ns_1) = -(n s -l) di.di.C~ ^ = -C-r 2+ ( n °~V 



6ij + (n s -3)-jj 



(4.45) 



We see that the result does not vanish for n s = 1. What vanishes for n s = 1 is the trace of this 
tensor. This means that for a scale invariant spectrum, the isotropic rescaling due to lensing does 
not contribute to the 3-point function. This makes sense in light of the discussion above: for the 
isotropic part there is no angular weight so that everything works as for the consistency relation for 
the squeezed limit of the primordial 3-point function [27, 41, 42]. On the other hand, the anisotropic 
case is similar to what happens when one calculates the primordial 3-point function of a graviton and 
two scalar modes, in the limit when the graviton wavelength becomes very long. The gravitational 
wave induces an anisotropic rescaling of the scalar 2-point function and the result does not vanish 
for a scale invariant spectrum [27]. An analogous effect is found when computing the contribution 
to the scalar trispectrum from graviton exchange [43]. In the limit where the graviton wavelength is 
very long, the non-isotropic rescaling induces a correlation between a pair of scalar 2-point functions. 
This effect has the same spin-2 angular dependence as the lensing. 

Let us go back to eq. (4.44). In our case the normalization of the spectrum is given by C = 
A/{9tt), so that the expression of the 3-point function in the squeezed limit gives 

{a f a t2 a ti ) = {2n) 2 S{h + f 2 + f 3 ) • *^<^>'(Z L ) (- A) -1 ^ _ 2^) . (4.46) 

The correlation between the temperature and the deflection angle is given by 

l — a .)'(l T ) - -11- [ — ^ Tdr T ~ T ^ A e ^(r-r.) - _L . AA 

(4.47) 



Thus we have 



B 



lens 



_ 4A 2 l Li l Lj f JsilsA U4R) 
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In the limit I2 — > 0, the explicit expression (4.35) gives, taking into account the permutation l\ «-> I3 



B 



lens 



AA 2 1 
27tt 2 Z| 



h -h h- (h + h) 



% {h + kf 



AA 2 l 2i l 2j d I 



27tt 2 i\ dhj q ' 



(4.49) 



which coincides with the expression above. 



5 The total CMB bispectrum 

In the previous section we have separated the calculation of the CMB bispectrum generated in the 
Sachs- Wolfe limit into five contributions: an intrinsic contribution expressed in terms of the New- 
tonian potential evaluated at last scattering, in eq. (4.6), the Rees-Sciama effect, in eq. (4.13), a 
contribution from the time dependence of the vector and tensor components of the metric, respec- 
tively in eqs. (4.20) and (4.28), and finally the lensing effect, in eq. (4.34). However, it is important 
to stress that only the sum of these contributions has a physical, gauge invariant, meaning. In this 
section we turn to discuss this sum, i.e. the total bispectrum. This is plotted in figure 8. 




0.2 0.4 0.6 0.8 1.0 

Figure 8: The total CMB bispectrum. 



By comparing this with figure 7 one can appreciate that the lensing effect largely dominates the 
total bispectrum 8 . Let us see this more quantitatively. 

8 As already stressed, the separation among different effects is gauge-dependent. Here we use the standard 
terminology in calling "lensing" the deflection of photons in Newtonian gauge. For a discussion about the 
gauge-dependence of lensing, see [39] . 
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In the squeezed limit the bispectrum is dominated by the intrinsic contribution and the lensing. 
In this limit we can compare the total bispectrum to the local bispectrum (3.16) taken with f^ al = 1- 
This yields 

TTi TT-, 7. TV = _ 1/6 — COs(2#) , (5.1) 

£ local (l,r 2 -> 0,r 3 -> 1) / v /> v y 

where represents the angle between the short and long wavelength modes l\ and / 2 - Thus, the 
total bispectrum corresponds to f^i? 1 = —1/6 — cos(2#). 

Note that this result can be obtained by simple arguments. As explained in section 2, the factor 
—1/6 can be inferred using the fact that a mode still out of the Hubble radius today cannot affect 
a physical measurement. The angular dependent factor — cos(2#) can be inferred by looking at the 
effect of a long wavelength lensing mode on the power spectrum, as explained in section 4.5. 

A remark on the angular dependence in eq. (5.1) is in order here. Although it is non-vanishing 
in the squeezed limit, the lensing contribution (4.34) is not of the local form (3.16). In particular, 
as the angular dependence averages to zero, a non-Gaussianity test based on a local estimator of 
the form (3.16) would be almost blind to the lensing signal. A quantitative way to measure how a 
signal overlaps with another is provided by the cosine between two bispectra, defined as [23] 

«»(Si,B 2 ) = rR^W B /R—FT ' (5 ' 2) 
where B\ ■ Bi is the scalar product between two bispectra, given by 

BiB2 = ^yi^p 6c^q 3 — (5 - 3) 

oc /dr 2 dr 3 ^rMWM^ (5.4) 

Indeed, we find that the cosine between the lensing bispectrum (4.34) and the local bispectrum 
(3.16) is cos(i?i cns , -Biocal) = 0.03. 9 For instance, one can compare this to the cosine between the 
local and equilateral bispectra, which is much larger, cos(-B e quii, -Biocal) = 0.30. Thus, due to the 
angular dependence of the squeezed limit, the lensing signal is orthogonal to the local one. We can 
now compare the total bispectrum to the local one. The cosine is cos(.B tota i, Biocal) = —0.17. Thus, 
as it is dominated by lensing, the total bispectrum is almost orthogonal to the local signal. However, 
due to the term — 1/6 in eq. (5.1) the orthogonality is not complete and the total bispectrum slightly 
overlaps with the local one. 

In the equilateral limit all the five contributions to the total bispectrum become important. 
However, the lensing numerically dominates. In this limit we can compare the total bispectrum to 
the equilateral bispectrum (3.19) taken with = 1. This yields 

B total (l, 1, 1 



B e< i uil (l, 1, 1) 



3.13 . (5.5) 



9 Note that the scalar product with the local bispectrum is logarithmically divergent for r-i — ► 1 or r% — ► 1. 
Thus, in order to evaluate it we have put the cutoff r max = 0.999. 
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Shape: 


total 


local 


equil 


lens 


total 


1.00 


-0.17 


0.41 


0.98 


local 




1.00 


0.30 


0.03 


equil 






1.00 


0.47 


lens 








1.00 



Table 1: Cosines between different shapes of bispectra. 

Thus, the total bispectrum corresponds to /ni!" 1 = 3.13. As it is not vanishing in the squeezed limit, 
its cosine with the equilateral shape will be smaller than unity. Indeed we find cos(-Btotai) -Bcquil) = 
0.41. Note that this value is larger than the cosine between local and equilateral shapes, i.e. 0.30. 
Thus, the total bispectrum is "more equilateral" than the local one. Finally, to have a confirmation 
that the lensing effect dominates the total bispectrum, we can compute the cosine between the total 
signal and the lensing. This is cos(B tota i, -Biens) = 0.98, which is very close to one, as expected. A 
summary of the cosines is given in table 1. 

It is important to stress that the shape associated with lensing, with an angle dependent squeezed 
limit, represents another interesting template for the bispectrum besides the local, the equilateral 
and the ones studied in [44, 45]. As it is rather orthogonal to the standard local and equilateral 
templates, in the future it would be interesting to put limits on it, even independently of lensing. 

6 Conclusion 

In this paper we have calculated, assuming perfect matter dominance, the complete CMB bispectrum 
on large angular scales, larger than the Hubble radius at recombination, considering for the first time 
all the relevant effects. Although our results give the exact bispectrum in a well defined physical 
limit, there are many ways to improve our calculations to make them closer to the real universe. One 
should include the recent dark energy domination and the early transition from radiation to matter 
dominance along the lines of [16]. This will give qualitative new phenomena, like the rather large 
ISW- lensing correlation [46]. Given that we are on large angular scales, a full-sky treatment would 
be more precise than our flat-sky expressions, although the results for the bispectrum will be much 
more complicated and difficult to understand. Finally, the small deviation from a scale invariant 
spectrum should be included. Taking all this into account would give the correct prediction for our 
universe of the large angle bispectrum. This is clearly far from the complete answer. The modes 
on scales larger than the horizon at recombination are quite few and most of the bispectrum signal 
comes from triangles with modes on sub-Hubble scales. Entering in a sub-Hubble regime requires 
the whole machinery of second-order Boltzmann equations that we have not touched in this paper. 
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The calculated bispectrum is rather small: the final bispectrum is dominated by the lensing 
contribution, which gives f^{] al = — cos(2#), with 9 the angle between long and short modes. Even 
if we could use our results on arbitrarily short scales, this would be below Planck sensitivity, limited 
to f^ al ~ 5. This means that the bispectrum in the Sachs- Wolfe limit does not represent a relevant 
contamination for the forthcoming searches for primordial non-Gaussianities. 

A way to go beyond the large angle regime is to correct the results of [22] to get the full 
bispectrum in the squeezed limit, with one (but not necessarily all) of the modes on scales larger 
than the horizon at recombination. We leave all these directions for future work. 
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Appendix 

A Flat-sky and integrated effects 

At first order, the gravitational contribution to the temperature anisotropies in matter domination 
is the Sachs- Wolfe effect, 



where n is the unit vector specifying the line of sight direction, D e = t q — r e is the (background) 
conformal distance to the last scattering surface and <p is the first order Newtonian potential. In the 
flat-sky formalism [31, 32], one chooses a fiducial direction z and expands at the lowest order in the 
angle 6 between z and n: 




(A.l) 



n = (sin 9 cos </>, sin 9 sin 0, cos 9) ~ (m x ,m y , 1) , 



(A.2) 
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to being a 2-dimensional vector normal to z. The multipole is simply the 2-dimensional Fourier 
transform with respect to to: 

a r = J d 2 m e-^^n) = \ J ^«5(f- k\\D e )e ik - D ^ % . (A.3) 

One can show that the flat-sky multipole corresponds to the large / limit of the full-sky one. The 
two are related by 



a l 



l7T -Y^r m a lm e m ^ , (A.4) 



21 + 

m 



'2Z + l. m 
ai m = \j — A 1 



47T 

Similar expressions hold also for the power spectrum and the bispectrum. The power spectrum is de- 
fined as (oj- a^) = (2tt) 2 8(11+12)0^ in flat-sky approximation, and as (ai imi a[ 2m2 ) = S^^S^Cf^ 
in full sky; the two expressions are related by Cf u11 ps Cf &t for large I. The bispectrum in the full 
and flat sky are defined respectively as 

(a ri a t2 a t:i ) = (2n) 2 6(h + h + h)B(h,T 2 , h) , (A.6) 

(ai )=\ B hhl3 , (A.7) 

\mi to 2 to 3 J 

where ( m \ }% 2 Jfi 3 ) is the Wigner 3-j symbol. The two expressions are related by: 



Bhhh a g y ^ (ai + 1) (^ +1)(a3 + 1 ) B(ri 4r3) (A 8) 

The derivation of these expressions can be found in [32]. 

To better understand what happens when we correlate effects which are important at different 
times, we can do a simple exercise 10 : we calculate the 2-point function of two integrated effects 
which peak at different times t\ and t 2 . We will see that the correlation decays exponentially when 
r ^ ( r 2 — r i)A> an d that the power spectrum is proportional to S(h + h) up to exponentially small 
terms. Consider a generic integrated effect at first order: 

a r = jT dr J ^<5(f- fc||(r - r))e ik ^~^ 9' (r)^ , (A.9) 

where g(r) is a growth function. Now we correlate two such effects, with different growth functions 



10 We thank F. Bernardeau for suggesting this example. 
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g(r) and /(r): 



(alaZ) = 4tt I dr a f{r a ) I dng'in) I d 3 kd(h - k\\{ To - T a ))5{l 2 + k\\{r - r 6 ))e 



ik ± (T b -T a ) A 



4vr / dr a 



f'(Ta 



dng'in) / dk±6 [h + h + h 



x e ik±(T b -T a ) 



A 



(kl + ll/(T -T a ) 2 ) 3 / 2 



For simplicity, we approximate the growth functions with step functions, such that 

/'(r) ~ S(t - ti) , g'(T) ~ S(t - r 2 ) , 



(A.10) 



(A.ll) 



where we consider T e < ti < t 2 < t c . Thus we find 



(44) 



r<5 Ni + «2 + <i 



dA; ± e^ (T2 - ri) (A:i + Z?/(t - n) 2 y 3/2 . (A.12) 



ii '2 (r - ti) 2 \ r - Tx 

The integration over /cj_ can be done analytically, yielding 

(aiaf) = (27r) d <i + / 2 + Z-i -7 r r ifi Z r 

!l i2 V T -TiJ 7T (T - Tl) Zl 



(A.13) 



t - n 

where i^i is the modified Bessel function, with asymptotic behaviours K\(x) — > l/x for x -C \/2 



and Ki(x) — > -\/tt/2x e~ x for s > 3/4. 
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Figure 9: The function xK\(x) with the x axis in logarithmic scale. 
We can see that, when Zi ^ 2 ~^ •> 1, the correlation is exponentially suppressed. Then, in the 



limit Zi <C (t — Ti)/(r2 — ti) we find 



h h ir If 



(A.14) 



In general, the translational invariance in 2d is only approximate; however, the approximation is very 
good since for large multipoles the correlations are exponentially suppressed if the sum ^ i Zj 7^ 0. 
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B Detailed calculation of the Rees-Sciama effect 



In this appendix we compute the Rees-Sciama bispectrum. We start from eq. (4.13) before summing 
over cyclic permutations: 

A 2 1 



B 



RS 



378vr 2 If 



I foo p-\-oo 

\ dxx d yi dy 2 e^+y^iyl + r\)-^\y\ + r 2 2 )~ 3 / 2 

1 J J -oo (B 1) 

3 ^i 2 + \r\ + 4 + 2y m + l(y\ + y 2 ) - 



2x 5 (yl -yl + r l- r lf 



.2 1 ' 2' ' J ' ™ ' 2^ ' 2 ( yi+2/2 )2 + r | 

Since there are some pieces in the kernel that can be integrated analytically, we now compute them 
as a check of our numerical integration. To proceed, we make use of the following known integrals: 



/+oo 
-oo 



dy- 



iyx 



-oo (y 2 + a 2 )a 

/+oo 
dye iyx 
-oo 

/+oo 
dye iyx 
-oo 



3- - — Ki(ax) , 



y 


. d 


(y 2 + a 2 )2 


dx 


y 2 


. d 


(y 2 + a 2 )2 


dx 



= 2ixKo(ax) , 
(2ixKo(ax)) = 2 [Ko(ax) — axK±(ax)] , 



(B.2) 
(B.3) 
(B.4) 



where the K{ are Bessel modified functions. We then split the bispectrum into four pieces, three of 
which are integrated analytically: 

2A 2 1 /3 



= - 



1897T 2 if 
£>( 2 ) 



L /3 3 \ 1 f°° 

j y 2 r i + 2 r2 + r v ^ y dx x3Rl ( riX ^ Ki ( r2X ) ' 



4,4 2 1 



189vr 2 if 



1 ■/ o 



dxx 3 K (na;)^o(^2^) , 



(B.5) 
(B.6) 



S (3) 



5A 2 1 f 00 



1897T 2 If 



1 •/ 



dx x 



1 1 

x ( — ^ri(rix)Ko(r2x) H K$(rix)Ki(r 2 x) 

r\ T2 



B^ = 



5A 2 1 f°° 



4 / dxx / 

1 JO J— oo 



756tt 2 If 

The first three pieces can be integrated in x, giving: 



-x 2 [ — + - ) K 1 (nx)Ki(r 2 x) 
n r2 



d yi d y2 e^+^( yi 2 + r 2 )-|( y 2 +r 2 )-| (y? 2/2 + r? ^ 



B (i) 



1 



44 2 1 /3 2 3 2 2 
1897 2 Zf V2 ri + 2^ + r V r 2 r 2 (r 2 - r 2 f [ 



(yi +y2) 2 + rl 



r\ — r\ — Ar\r\ In — 



16A 2 1 



189vr 2 if (r 2 - r 2 f . 



r\-r\- {r\ + r 2 ) In 



ri 
r2 



5A 2 1 



189vr 2 if (r 2 - r|) 3 



4 4 

To r 



5( r 2 - r 2 ) + 1| _ 11 + 2 (r 2 - r 2 2 ) In ^ - 2(3r 2 + 5r 2 ) In ^ 



T2 



T2 



n 



(B.7) 
. (B.8) 

(B.9) 

(B.10) 
(B.11) 
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The fourth piece (B.8) cannot be integrated analytically, making the numerical integration necessary. 
However, comparison between the numerical integration of the other three pieces and the analytical 
expressions (B.9) to (B.ll) gives consistent results. This provides a check of the validity of our 
computation. 

Now we turn to the numerical integration of eq. (B.l). To do it we first have to perform 
analytically the x integral, which is ill-defined. In order to overcome this problem we first change 
variables from y±, y2 to y + = y\ + y 2 and y_ = y\ — y 2 , and then regularize the integral in the 
following way: 

f + OO r + OO f+OO /"OO Q 



poo r+oc r+oo r+oo r+oo poo 

/ dx dy + / dy-Xe %y+x f(y + ,y-)=-i / dy + / dy_ / dxf(y + ,y 

JO J —oo J —oo J —oo J —oo JO 



ty+x 



dy^ 

/+oo r+oo Q roo 

dy + / dy_— f(y+,y-) / dxe iy+x 
-oo J-oo dy+ J 

/+oo r+oo i a 

dy+ / dy_ — — f(y + ,y-). (B.12) 
oo J-oo V+ dy+ 

The last integral follows from the prescription 

dxe iy+x e- £X = . (B.13) 

-iy+ + £ 

After integrating in x, we obtain: 

A 2 1 f +QO , , Id' 



B RS = 



378vr 2 I 



i / °° ^V2-^~ \{yl + rl)-*/\yl + rl)-^ 
,f J-oo y+ dy+ 



3 2 .322, 2 „ 17J _ 1 ^„2 , 7 .2\ _ 5 (3/1 ^2+ r l r l) 21 



(B.14) 



which, after changing variables from (2/1,2/2) to (j/ + , y_), and performing the derivative, gives a form 
which can be integrated numerically. The final results of the integration, after summing over cyclic 
permutations and setting n = 1, are presented in figure 4. 

We can compare the Rees-Sciama contribution with the local shape in the squeezed limit. We 
will see that while the local shape diverges as 1/r 2 in this limit, the Rees-Sciama only diverges as 
1/r. Going back to eq. (B.l), we can study the behavior when one of the r goes to zero; notice that 
the expression must be symmetrized so that we have to study both the limits r 2 — > and r% — > 
in eq. (B.l). For r 2 — ► we have an infrared divergence in the y 2 integral coming from the power 
spectrum which goes as y^ 3 for r 2 = 0. This would give a divergence r^ 2 as in the local model. 
However, for r 2 = and r\ = = 1 the expression in brackets in the second line of (B.l) goes as y 2 
for j/2 — ^ 0, but its integral vanishes due to parity, leaving only terms which are at most logarithmic 
divergent and thus suppressed with respect to the local shape. An additional divergence comes from 
the limit — > in eq. (B.l); in this case the integral diverges in the limit yi + y 2 — > 0. Notice that 
in this case one also has to take into account the integral over x which diverges for yi + y 2 = 0. To 
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study the behavior for r% — > one must first integrate in x using the prescription (B.12). One can 
see that the leading divergence in the resulting expression comes from a term of the form r%/y+, 
which gives a I/V3 divergence. This is dominant compared to the divergence in r2, but it is still 
subdominant compared to the local case. We conclude that the Rees-Sciama result is subdominant 
compared to the local shape in the squeezed limit. This analysis is a good check of the numerics, 
which indeed shows a 1/r divergence in the squeezed limit. 
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